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Abstract 

We analyze the effect of polydispersity in the arm number on the effective in- 
teractions, structural correlations and the phase behavior of star polymers in 
a good solvent. The effective interaction potential between two star polymers 
with different arm numbers is derived using scaling theory. The resulting ex- 
pression is tested against monomer-resolved molecular dynamics simulations. 
We find that the theoretical pair potential is in agreement with the simulation 
data in a much wider polydispersity range than other proposed potentials. We 
then use this pair potential as an input in a many-body theory to investigate 
polydispersity effects on the structural correlations and the phase diagram 
of dense star polymer solutions. In particular we find that a polydispersity 
of 10%, which is typical in experimental samples, does not significantly alter 
previous findings for the phase diagram of monodisperse solutions. 
PACS 82.70.Dd, 64.60.Fr, 61.20.Ja, 61.41.+e 

I. INTRODUCTION 

A star polymer consists of linear polymer chains tethered to one common central core. 
The number of polymer chains / is usually referred to as arm number or functionality. If 
the degree of polymerization, i.e., the number of monomers per chain, N, is the same for all 
chains, the star polymer is called monodisperse with respect to N or "regular" . For large N, 
the size of the central core particle is much smaller than the overall extension of the whole 
star. For a single star, the density of monomers around the central particle, </>(r), is radially 
symmetric, and from this density profile the so-called corona diameter a, which measures 
the extension of the star, can be defined as the diameter of a sphere around the star center 
where all Nf monomers are found "inside" IJ. 

In the last years, star-shaped polymer aggregates have attracted a considerable amount 
of interest from both the experimental and the theoretical point of view |2|-|J- The reason 
for that is threefold. First, from a technical point of view, star polymers are important for 
several industrial applications [0. One example are hydrogenated polyisoprene star poly- 
mers, which are used as viscosity index modifiers in oil industry applications due to their 
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excellent shear stability. Further, commercial star polymers are brought into action in coat- 
ing materials, as binders in toners for copying machines, and in several pharmaceutical and 
medical applications 0. Second, from an experimental point of view, the recent synthesis of 
regular star polymers with various possible numbers of arms by Roovers and coworkers |5]|| 
made it possible to explore the physics of well-defined model systems, which are monodis- 
perse in both the number of arms and the degree of polymerization. Important examples 
are polyisoprene stars with / = 8, 18 M and polybutadiene stars with / = 32,64, 128 ||, 
both synthesized by anionic polymerization. Third, star polymers constitute an important 
soft-condensed matter system, linking the fields of polymer physics and colloid physics, thus 
attracting also interest from a purely theoretical point of view. Star polymers with small arm 
numbers (/ = 1, 2) resemble linear polymers. Thus, their configurations show a consider- 
able asphericity [|S|-|iO|j, although their chain-averaged number density of monomers, 0(r), is 
spherically symmetric around the center of mass of the polymer. With increasing arm num- 
ber /, the asphericity of the stars has been shown to decrease considerably |TT|-|ISJ, leading 
to "stiff" spherical particles in the limit of high /. It is essentially this limit that a descrip- 
tion of star polymers as sterically stabilized colloidal particles holds. This polymer-colloid 
hybrid character of star polymers has been explored in a number of publications dealing 
with the structural P,[7|j2l]|-|52]] and dynamical [|2l5|-f2l|P3}-|3"5[ properties of star polymers. 

In the usual theoretical description, star polymers in dense solution are considered as an 
effective one-component (monodisperse) system, i.e., all the stars have the same arm number 
/ and all the linear polymer chains attached to the center have the same molecular weight. 
While the latter can be realized by a careful preparation method, the former is in general not 
true in actual samples. The preparation of star polymers proceeds by adding linear polymer 
chains to a dendrimer-like core with reactions centers placed at the end of the dendrimer 
and the polymer chains. The chemical reaction is typically incomplete such that few of 
the reaction centers at the dendrimer are not linked to a linear chain. Hence an intrinsic 
polydispersity in the arm number arises. By ultracentrifugation one can estimate the relative 
polydispersity in the arm number to be around 5-15 percent in carefully and slowly prepared 
samples. It can, however, also be much larger for fast reactions |3^,f|(J. In a comparison 
between experimental data and theory, in almost all previous studies the assumptions of a 
monodisperse sample was implicitly made such that a natural question concerns the influence 
of polydispersity on the statistical properties of star polymer solutions. As far as theory is 
concerned, the situation resembles much the case of colloidal suspensions where the effect of 
charge and size polydispersity has been a topic of intense recent research, see Refs. [CT-flH 
for recent reviews. 

In this paper, we investigate the effect of arm number polydispersity on the effective 
interaction, the structural correlations and the phase diagram for star polymer solutions in 
a good solvent. Our work is based on a theoretical analysis using scaling theory and computer 
simulations. We first derive the effective interaction between two star polymers of different 
arm numbers f\ and from scaling assumptions. Basically the interaction is logarithmic 
with the core-core distance r between the two stars but the prefactor depends explicitly on fi 
and /2 and differs from that of a monodisperse sample. The resulting pair potential is tested 
against molecular simulations and good agreement is found even when f\ is very different 
from /2 while earlier descriptions |^7j are found to hold for small polydispersities but to be 
too simplistic for large polydispersities. In a further step towards a full description of dense 
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star polymer solutions, we then perform computer simulations of a classical many-body 
systems interacting by means of this effective pair interaction. We find that a polydispersity 
in the arm number reduces the structural correlations but this effect is less drastic than 
for size polydispersity in sterically stabilized colloidal suspensions. Further, we present 
evidence that the reentrant melting behavior predicted from the theoretical treatment of a 
monodisperse description [pl], [3^ , [31] , |48|l does not change drastically in a polydisperse solution 
of star polymers. 

The paper is organized as follows: In section |TJ we present the scaling ideas and derive 
the effective interaction potential between two stars of different arm numbers. A test of this 
result against molecular computer simulation is performed in section [TTT[ Results for the 
structure and phase diagram from a simulation of a many-body system are given in section 
prV] . Finally, we conclude in section |V|. 



II. SCALING THEORY FOR AN EFFECTIVE POTENTIAL 

The effective interaction between two star polymers results from the steric interaction 
of all monomers that constitute the polymer chains on each of the stars. It is obviously a 
formidable task to derive the effective potential from first principles. Luckily we are in the 
position to give a rigorous result for the limiting case of short distance between the two star 
cores, while good arguments exist for the description of the long range part of the interaction. 
We have found in the treatment of the interaction of monodisperse star polymers that the 
combination of these two approaches leads to good agreement with simulation results as well 
as with scattering experiments on real star polymer systems 



We first derive the short distance interaction using scaling arguments. Many details of 
the behavior of polymer solutions may be derived using the renormalization group (RG) 
analysis [[51]. Here, we use only the more basic results of power law scaling: the radius of 
gyration Ri(N) of a single linear polymer chain with N monomers and its partition function 
Zi(N) are found to obey the power laws: 

Ri(N) ~ N v and Z X {N) ~ z N N~ vri2 . (1) 

The fugacity z measures the mean number of possibilities to add one monomer to the 
chain. It is microscopic in nature and will depend on the details of the theoretical model 
or the experimental system. The two exponents v and r]2 on the contrary are universal 
to all polymer systems in a good solvent, i.e., excluding high concentration of polymers or 
systems in which the polymers are collapsed or are near the collapse transition. A family 
of additional exponents r]f governs the scaling of the partition functions Zf(N) of polymer 
stars of / chains each with N monomers while the power law for the radius of gyration 
remains unchanged: 

R f (N) ~ N v and Z f (N) ~ z fN N v ^- fr,2 \ (2) 

The exponents of any other power law for more general polymer networks are given by 
scaling relations in terms of rjf and v [p^-|o^]. Note that eqs. (|I]) and (FJ) imply rji = 0. The 



values of these exponents are known from RG and Monte Carlo (MC) simulations |55|. In 
the limiting case of many arms, / ^> 1, the leading behavior of r)f is given by j56 
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Vf - -f>\ 



(3) 



a result that has been found also from of geometrical considerations using the blob model 
of a star polymer pD |. The scaling law for the partition sum of two star polymers may 



be derived from a short distance expansion 57,58]. The partition sum of the two stars 



2/i/ a (iV", r) at small center-to-center distances r factorizes into a function Cf x f 2 (r) and the 
partition function Zf 1+ f 2 (N) of the star with fx + f 2 arms that is formed when the cores of 
the two stars coincide: 

Z flh (N,r) ~ C flf2 (r)Z fl+f2 (N) . (4) 

Such a behavior is generally assumed in scaling arguments . At large star-star separation 
Zfaf^N, r) is the product of the single star partition functions. Again we may argue that the 
ratio of the partition function of two stars at finite separation r to that at infinite separation 
can only be a function B of r/R\ noting that R\ is the only relevant length scale in this 
problem {Rf x and Rf 2 differ from R\ only in an /-dependent prefactor): 

Z hf2 (N,r)^B flh (r/R 1 )Z h Z h . (5) 

Taking this into account and inserting the power law scaling according to Eqs. ([I]) and (0) 
we find that Eqs. (f|) and (||) are only compatible if also C{r) follows a power law 

C hh (r)~r e f^, (6) 
and the so-called contact exponent QVt 2 obeys the scaling relation 

Q/ik = Vfl + Vh ~ Vfl+h ■ ( 7 ) 

Note that in terms of the RG of polymer field theory the exponents T]f correspond to dimen- 
sions of operators associated with the partition functions of single stars [|54]] and relation 



(0) is a direct consequence of the short distance expansion. The mean force Ff x f 2 (f) acting 
on the centers of two star polymers with f\ and f 2 arms is now derived as the gradient of 
the effective potential V eS (r) = —kBTlog[Zf 1 f 2 (r)/(Zf 1 Zf 2 )) with k-^T denoting the thermal 
energy. For the force at short distances r this evidently results in 

(s) 

F flf2 (r) = k B T^. (8) 
Using the above many arm limit for r/j one may match the contact exponents OVi to the 



known values for / = 1, 2 [|59|] fixing the otherwise unknown prefactor in eq. (H). Assuming 



(sT 

that the behavior of the may be described by this approximation for all / one finds: 

Fff(r)*k B T- J — . (9) 

This matching in turn suggests an approximate value for the rjf exponents, taking into 
account rji = and r] 2 ~ —5/18: 
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Inserting this into the general formula for the interaction of two different star polymers the 
contact exponent reads for large f\ and ^ 
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(A 3/2 + fT) 



(11) 



We note that on a phenomenological level there are other possibilities to describe the 
interaction in a polydisperse system of star polymers. In colloidal systems two other ap- 
proaches to the interaction of polydisperse particles are commonly used: in a steric stabilized 
system with a polydispersity in the radii of the particles one expects to find an interaction 
radius that is the arithmetic mean of the radii of the two interacting particles [4*4",|B"0"|. In a 
charge stabilized solution the other hand, the interaction of two particles is proportional to 
the product of the two charges and the effective charge is calculated as the geometric mean of 
the two charges. The latter description was used in an earlier investigation of star polymers 
treating / as an "effective charge" parameter p7| , |6T1 , |62 
possible approaches we define: 



For comparison with these other 



7l + /2 
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(12) 
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We now turn to the interaction at larger separation. It has been recently shown that 
a Yukawa like tail for the potential between star polymers reproduces the results of both 



simulation and scattering experiments for monodisperse solutions |28 ^9| , [50[| . The natural 
scale for a star polymer is the corona diameter o~f = 2XRf with A ~ 2/3 |49|]. For the 
potential between two /-arm star polymers at distance r > o~f one has 

V ff {r) ~ -exp(-r«/). (14) 

The decay length I/kj = 2a / / v7 is the diameter of the outermost blobs of the star polymer 
in the Daoud-Cotton model jl|j28|. 

In a simple but successful approach the full potential for the interaction between two 
/-arm star polymers was constructed by concatenating the long range and short range 
potentials. In this model the crossover from short to long range behavior takes place when 
the distance of the two cores is o~f. It is natural to generalize this to two star polymers with 
/i and arms in such a way that the crossover point is at distance a — (o^ + <7/ 2 )/2. It 
is less obvious how to choose the decay length of the long range part. Our choice 1/k — 
l/Kfr + l/Kfa is in accordance to the simulation results presented below. Note that although 
we have omitted the subscripts on a and k, the latter are functions of /i and /2. 

The full potential then reads 
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T 1 

ln(— ) H for r < a; 



, ir -VftA(r) = ejf A !»„ l + « (15) 

B expfcTK — r/c) tor r > a. 

I 1 + oTcr v y 

In the monodisperse case /x = f 2 this potential reduces to the one between two identical 
star polymers that has been successfully tested in extensive simulations. 



III. COMPARISON WITH SIMULATION RESULTS 

We have performed Molecular Dynamics (MD) computer simulations of two star poly- 
mers with different arm numbers f\ and f 2 by resolving the monomers as classical beads 
along the chains. We use the same simulation model as in previous studies of star polymer 



solutions in a good solvent [£5J63"f. The arm length, i.e., the number N of beads along each 
chain, was kept fixed to iV = 50, 100 which was shown in previous work to be large enough 
to guarantee a sufficient scaling behavior. The main features of the simulation model can 
be summarized as follows: 

1. A purely repulsive truncated Lennard- Jones interaction between all iV(/i + f 2 ) 
monomers is assumed. 

2. An attractive FENE-potential is added for the interaction between neighboring mo- 
nomers along a chain. 

3. To accommodate the polymer arms, a hard core with radius is introduced at the 
center of the stars. The cores interact with the monomers according to the above 
mentioned potentials with a separation shift of R^ d \ 

Two stars with arm numbers f\ and f 2 were fixed at core separation r. After a long equi- 
libration period, the total force F^f 2 acting onto the cores is averaged. The temperature 
T is kept constant. Both the equilibration time and the time during which averages were 



performed was about lOOOr where r = ymcr 2 /e is the Lennard- Jones time unit with a, e 
and m denoting the Lennard- Jones length, energy scale and the monomer mass. A simula- 
tion snapshot is shown in Fig. |]. When comparing the mean force to the derivative of the 
theoretical effective potential, two technical problems arise: (i) In contrast to the theory, we 
have a finite core size R^ d > which is the same for both stars, (ii) The two corona diameters, 
°7) (/ = iL/2) which enter in the potential of eq. fllS]), are not directly accessible in a 
simulation. 



As in the monodisperse case the first difficulty (i) is overcome by plotting the 
inverse force 1/Ff 1 f 2 versus distance r. The divergence in the force occurs already at core 
separations r = 2R^ which leads to a zero in the 1/Ff 1 f 2 plot. The inverse slope of the 
l/F^-plot then corresponds to the prefactor 6j*j 2 of the logarithmic term in the effective 
pair potential. It is this slope which can directly be compared to our theoretical prediction 
eq. (|T5|). Furthermore, a linear function in the 1/Ff l f 2 plot is a direct check for the validity 
of the In r-term in the effective interaction potential, eq. ([15]) . 

To handle the second difficulty (ii) we consider the radius of gyration Rf of each star 



which is readily accessible in a simulation of a single star ||49|| . We assume a proportionality 
between the Rf and <j/, 
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a f = 2XR f 



(16) 



and use A as a fit parameter to fit the full force versus distance curve for given arm numbers 
fx and fi. The results for the optimal fit parameter are shown in Table |. We obtain an 
averaged value of (A) ~ 0.66, which is nearly independent of fx, / 2 . Thus, the results found 
here are consistent with those from the previously investigated monodisperse case |49(| . We 



further remark that this value coincides with that used in Ref. |28| to fit experimental data 
for / = 18 and that A is independent of N, consistent with scaling theory 

Results for the comparison between theory and simulation are shown in Figure where 
the reduced inverse mean force l/F^fa between the star centers is plotted versus intercore 
distance r for two arm numbers fx = 10 and f 2 = 50. We observe two important facts: 
first, the data indeed fall on a straight line proving the logarithmic behavior of the potential 
inside the mean corona. The straight line hits the origin at r = 2R^ showing the relevance 
of the finite core in the simulations. Second, the inverse slope Qf ± f 2 agrees very well with 
the theoretical prediction also shown in Figure |2|. A similar behavior was observed for all 
other combinations of arm numbers contained in Table |I[ 

In Fig. [| the prefactor O^f? is plotted versus / 2 , at a fixed value fx = 10, and compared 
to the different theoretical predictions in eqs. (p|)-([T3|). Of course, in the monodisperse 
case fx — f2 — 10, all the different theoretical expressions for coincide and agree 



well with the simulation, consistently with our earlier work [49]. For increasing asymmetry 



between fx and f 2 , significant deviations between the predictions of (a) an arithmetic or 

»(■) 

) 

'/1/2 



(g) a geometric mean and the simulation results become visible. While the expression 0^ 
gained from scaling theory is able to describe the slope even for large asymmetries, OV%„ is 



worse at large asymmetries while ©yw 2 even possesses the wrong curvature as a function of 
/ 2 . Similar conclusions hold for other combinations of fx and ,/ 2 . 

In Fig. [|, we show the force Ff x f 2 versus distance r — R^, scaled by Rx2 = (Rf x + Rf 2 )/2, 
comparing the theoretical force with simulation results for three cases: fx = 10, 18, 30; f 2 = 
50; N = 50. Note that formally the core size vanishes in the theory. There is a good 
overall agreement inside the corona region between theory and simulation even for large 
asymmetries of fx and / 2 . Results for distances r outside the corona, r > a, are presented 
in Figure [5]. Here, we test the exponential decay of the force Ff 1 f 2 with distance r by 
plotting the logarithm of the force versus distance. The crossover of the inner-core data to 
a straight line outside the core is clearly visible. Within the simulation error bars, the slope 
is consistent with the theoretical prediction eq. (|15|). 



IV. POLYDISPERSITY EFFECTS IN DENSE STAR POLYMER SOLUTIONS 

To explore the effects of arm number polydispersity on the structure and phase behavior 
of dense solutions of star polymers, we further performed large scale Monte Carlo simulations 
using the pair potential approach as given by the potential in eq. (|T5|). Thus, on the one 
hand, this approach closely corresponds to our recent work on the structural and phase 
behavior of monodisperse star polymer solutions |29^pT| , and, on the other side, is quite 



similar to investigations of polydispersity effects in solutions of hard spheres |4^j60|j6^j65 
and charged colloidal particles P5|J61]J6^J56| available in the literature. 
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In our MC simulations, representative particle configurations of approximately 2000 star 
polymers in a cubic simulation box were created in the following manner. First, a starting 
particle configuration was built up by placing the particles on random positions in the box, 
each particle assigned with an arm number / chosen from a Gaussian distribution g(f) of 



mean / and variance p = y f 2 — f 



9(f) = —j==exp 
pyzn 



2p 



1 / -\ 2 



(17) 



We use eq. ( PD above to define the polydispersity p in what follows. The Gaussian dis- 
tribution chosen to describe arm number polydispersity agrees well with experimentally- 
determined arm number distributions of star polymers synthesized by anionic polymerization 
SOfl . Note that nonphysical, negative arm numbers, which are in principle not prohibited by 



a Gaussian arm number distribution, did not occur in the simulations for the values of / and 
p chosen here. Second, an equilibration phase of approximately 50 000 MC cycles was per- 
formed by allowing both translational particle moves according to the standard Metropolis 
scheme, and particle exchanges of randomly chosen particle pairs, again using the Metropolis 
rule to decide an exchange to happen or not. For the calculation of pair potential energies, 
eq. (|15[ ) was used, together with the scaling relation [|TJ: 



a h \fo) 



(18) 



After the equilibration phase, approximately 100 000 MC cycles were simulated to gather 
statistics for both the radial center-to-center distribution function g(r) and the center-to- 
center structure factor S(k) pi| , p2]P^ , |io]jD7| . In what follows, we will use a = aj as the 
basic length scale, and the mean packing fraction fj = p<x 3 7r/6 (p being the number density 
of the stars) as a measure of the density. 

In order to investigate the polydispersity effects on the liquid structure of dense star 
polymer solutions, simulations with varying polydispersities p were performed at (average) 
arm numbers / and mean packing fractions fj, known to correspond to the liquid state in the 
monodisperse case p = [0,^l|]. We show typical results for g(r) and S(k) in fig. ||. As can 
be seen, increasing the polydispersity leads to decreasing spatial correlations in the fluid, 
indicated by a decreasing principal peak height of both g(r) and S(k). Thus, as expected, 
the effect of arm number polydispersity in star polymer solutions is quite similar to the effect 
of charge polydispersity in solutions of charged colloids fl5l]66fl . Note that the anomalous 



behavior of S(k) reported for monodisperse star polymer solutions in Ref. |29 also holds for 
polydisperse star polymer solution for all polydispersities p < 14 simulated. 

In order to explore the evolution of the freezing and reentrant melting phase transitions 



found for monodisperse star polymer solutions |3Cj],[JTJ with increasing polydispersity, we 
employed the following strategy. We performed a number of simulations at state points cor- 
responding to the solid phase in the monodisperse case, gradually increasing the parameter 
p. In particular, we chose / = 40 and fj = 0.5 in a first set of simulations. Our results 
for g(r) and S(k) are given in Fig. |7|. Let us begin with the discussion of the g(r)-data. 
As the monodisperse case p = corresponds to an solid state of bcc-symmetry [j30| , |31]1 , the 
shown g(r) is, in fact, the radially-averaged pair correlation function of the stars in the 
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crystal state.Q Again, increasing p leads to decreasing structural correlations between the 
particles indicated by a decreasing principal peak height and a "smearing-out" of the sub 
peak-structure of g(r) seen for p = 0. This scenario also manifests itself in the structure 
factors shown in fig. 0(b). Notice that the corresponding structure factor for p = exhibits 
strong Bragg peaks, again indicating the monodisperse system to be in the crystal state for 
f = f = 40, f] = 1] = 0.5, and is not depicted for that reason. For nonzero polydispersities 
the simulated structure factors do not show Bragg peaks and the value of the principal peak 
heights S'max decreases. Thus, increasing the polydispersity in a sample of star polymers at 
f — f — 40, fj — T) — 0.5 leads to a melting transition from a bcc crystal to a fluid phase. 
We therefore expect polydisperse star polymer solutions to show an enlarged fluid phase 
region as compared to the phase diagram of monodisperse star polymers in Refs. |30| , |31~ll . 
Furthermore, the stability range of the fluid phase is expected to increase with increasing the 
polydispersity p. To corroborate this expectation in more detail, we have performed further 
MC simulations of polydisperse star polymers for various state points (/, fj), all correspond- 
ing to state points close to the freezing or reentrant melting line in the monodisperse case 
P0|]jn|| . For the ranges of polydispersity examined here, (p < 14), no change in the topology 
of the phase diagram and no new crystalline phases were found as a result of polydisper- 
sity. A quantitative calculation of the full phase diagram by more sophisticated simulation 
methods |64],|65| is beyond the scope of this paper and will be left for future studies. 



V. CONCLUSIONS 

In conclusion, we have analyzed the effect of arm number polydispersity on the effective 
interaction and the structural correlations between star polymers in a good solvent. An 
analytical expression for an effective pair potential, given by eq. (|15D, was put forward and its 
validity was confirmed by Molecular Dynamics computer simulation. This pair potential was 
subsequently used to simulate the structural correlations between many stars. As expected, 
correlations decrease with increasing polydispersity. At the same time, though, the effect is 
much less pronounced than for size-polydisperse hard spheres. This can be seen in analogy to 
colloidal suspensions where size polydispersity of sterically-stabilized suspensions is known 
to lead to much more pronounced effects than charge-polydispersity in charged suspension. 
The microscopic reason for that is that the effective interaction is much softer for charged 
colloids which is similar to our case of star polymers. 

We finish with a couple of remarks: first it would be interesting to compare, both qual- 
itatively and quantitatively, our theoretical predictions to experimental data. In fact, the 
intrinsic polydispersity in the arm number can be measured and structural correlations are 
accessible, e.g., by neutron scattering of a star polymer solution with marked cores. Second, 
it would be interesting to map out, for a given polydispersity, the full phase diagram in- 
cluding freezing into different solid structures. Also, it would be interesting to develop and 



1 In fact, from the typical particle-particle distances read off from the pair correlation function, it 
can be concluded that the crystal structure is bcc, a result which is in agreement with calculations 
of bond order correlation functions [31]. 
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apply a liquid integral equation theory to predict structural correlations in a polydisperse 
star polymer solutions. The output of such a theory could be checked against our computer 
simulation data of section [IV]. Third, we have not considered a polydispersity in the length 
of the linear polymer chains attached to the centers. It would be interesting to study this 
theoretically and compare the results to samples prepared in such a way that different linear 
chains are brought to the reaction centers of the dendrimers. Work along these lines is left 
for future studies. 
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TABLES 



fl 


h 


N 


A 


3 


5 


100 


0.511 


5 


10 


100 


0.643 


5 


18 


100 


0.667 


10 


18 


50 


0.675 


10 


30 


50 


0.673 


10 


50 


50 


0.692 


18 


30 


50 


0.696 


18 


50 


50 


0.668 


30 


50 


50 


0.714 







TABLE I. List of the simulated /i-/2-combinations and the corresponding results for 
A = a/(Rj 1 + Rf 2 ). N is the associated monomer number per arm. 
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FIGURES 



FIG. 1. Snapshot of two interacting star polymers with f\ = 18 and /2 = 10 arms. 

FIG. 2. Reduced inverse force ksT '/ (Ri2Ff 1 f 2 ) between the centers of two star polymers (for 
fi = W,f2 = 50 and N = 50) versus reduced distance r/R±2, where R12 = (R^ + i?/ 2 )/2. The 
error bars were obtained by averaging over the results of 10 independent simulations. The dashed 
line is a linear regression of the data. The solid line is the theory from chapter |n[ 

FIG. 3. Simulation results for the prefactor Qf±f 2 for fi = 5 . . . 50 and f\ = 10 fixed, in 
comparison to three different theoretical predictions. 

FIG. 4. Simulation results (symbols) and theoretical results (lines) for the reduced effective 
force R^F^fe/ftsT versus reduced distance (r — 2R^) / R\2- 

FIG. 5. Logarithm of the reduced force, In (72i2-F/ 1 / 2 /fc_BT), versus reduced distance, 
(r - 2RW)/Rx2, for ft = 10, f 2 = 18 and N = 50. The error bars were obtained by averag- 
ing over the results of 10 independent simulations. 

FIG. 6. (a) Radial distribution function g(r) and (b) center-to-center structure factor S(k) 
for polydisperse star polymers in the fluid state. In the MC simulations, the average arm number 
has been chosen as / = 32, the mean packing fraction as fj = 0.5. Polydispersities p as indicated 
in the figure. 

FIG. 7. (a) Radial distribution function g(r) and (b) center-to-center structure factor S(k) 
for / = 40, f\ = 0.5. Polydispersities p as indicated in the figure. 
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von Ferber et al., Fig. 2 




von Ferber et al., Fig. 4 
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von Ferber et al., Fig. 5 




von Ferber et al., Fig. 7 b 
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